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I .  SUMMARY 

This  report  describes  the  operation,  maintenance  and  research  activities  at 
the  Norwegian  Seismic  Array  (NORSAR)  for  the  period  1  October  1981  to  81  March 
1982. 


The  uptime  of  the  NORSAR  online  detection  processor  system  has  averaged  91.9%, 
'which  represents  an  improvement  of  2%  compared  to  the  previous  period.  Almost 
half  of  the  downtime  was  due  to  problems  in  connection  with  the  termination 
of  the  NORSAR  ARPANET  circuit.  The  SPS  has  been  working  well  and  has  caused 
less  than  1%  downtime. 

A  total  of  1681  events  were  reported  in  this  period,  giving  a  daily  average 
of  9.2  events.  The  number  of  reported  events  per  month  varies  from  210  in 
February  to  363  in  December. 

There  have  been  no  major  breakdowns  on  the  communications  lines,  but  some 
lines  have  had  periods  with  bad  performance,  especially  the  02C  line. 

The  new  DP  system  to  be  implemented  on  IBM  4331/MODCOMP  is  now,  except  for 
some  EOC  functions,  ready.  The  communication  problem  between  IBM  4331  and 
Modcomp  has  been  solved,  and  test  runs  have  been  done  with  no  errors.  Some 
improvements  in  programs  and  some  more  testing  still  needs  to  be  done  before 
the  new  DP  is  ready  to  replace  the  current  system. 

The  research  activity  is  briefly  described  in  Section  VI.  The  subsection  1 
discusses  a  method  for  simultaneous  determination  of  source  parameters.  Sub¬ 
section  2  presents  a  processing  package  for  on-line  analysis  of  data  from 
small-aperture  arrays.  Subsection  3  is  a  study  on  optimization  of  seismic 
array  configuration.  Subsection  4  and  5  describe  the  upper  mantle  seismic 
heterogeneities  beneath  Fennoscandia  and  Eastern  Europe.  Subsection  6  describes 
a  work  for  mapping  the  crustal  structure  and  thickness  in  the  coastal  areas  of 
southern  Norway  and  the  general  features  of  the  Moho  near  and  in  the  Oslo 
Graben.  Subsection  8  is  a  noise  study  at  various  sites  in  Finland,  in  northern 
Norway  and  in  southeastern  Norway. 
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II.  OPERATION  OF  ALL  SYSTEMS 

II .1  Detection  Processor  (_DP)^  Opera 1 1  o n 

There  have  been  178  breaks  in  the  otherwise  continuous  operation  of  the 
NORSAR  online  system  within  the  current  6-month  reporting  interval.  The  SPS 
is  still  working  well  and  has  caused  only  0.8%  downtime.  A  software  problem 
that  occurred  in  conjunction  with  stopping  the  online  transmission  of  data 
via  the  ARPA  net  caused  almost  half  of  the  downtime  in  the  period.  The 
uptime  percentage  for  the  period  is  91.9  as  compared  to  89.0  for  the  previous 
period . 

Fig.  II. 1.1  and  the  accompanying  Table  II. 1.1  both  show  the  daily  DP  downtime 
for  the  days  between  1  October  1981  and  31  March  1982.  The  monthly  recording 
times  and  percentages  are  given  in  Table  II. 1.2. 


The  breaks  can  be  grouped  as  follows: 

a)  SPS  malfunction  20 

b)  Stops  related  to  possible  program  errors  74 

c)  Maintenance  stops  5 

d)  Power  jumps  and  breaks  5 

e)  Hardware  problems  61 

f)  Magnetic  t ipe  and  disk  drive  problems  11 

g)  Stops  related  to  system  operation  2 


The  total  downtime  For  the  period  was  353  hours  and  38  minutes.  The  mean- 
time-between-failures  (MTBF)  was  1.3  days  as  compared  with  0.9  days  for  the 
previous  period. 


J.  Torstveit 


Detection  Process 
31  March  1982. 
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Online  System  Performance 
1  October  1981  -  11  March  1981 


II. 2  Event  Processor  Operation 

In  Table  II. 2.1  some  monthly  statistics  of  the  Event  Processor  operation  are 
given : 


Teleseismic 

Core  Phases 

Sum 

Daily 

Oct 

81 

186 

13 

239 

7.7 

Nov 

81 

174 

19 

233 

7.8 

Dec 

81 

243 

120 

363 

11.7 

Jan 

82 

292 

58 

no 

11.3 

Feb 

82 

168 

42 

210 

7.1 

Mar 

82 

231 

15 

286 

9.2 

1294 

387 

1681 

NO 

NJ 

TABLE  1 1. 2.1 


B.  Kr.  Hokland 
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I 1.3  NORSAR  Data  Processing  Center  (NDPC)  Operation 
Data  Center 

The  operator  function  is  now  being  carried  out  by  four  persons  on  a  weekly 
shift  turn.  In  addition  a  fifth  person  steps  in  for  a  day  or  two  when  neces¬ 
sary,  as  they  all  have  another  main  function  to  cover. 

J.  Torstveit 


II .4  Array  Communication 

Table  II. 4.1  reflects  the  performance  of  the  communications  system  throughout 
the  reporting  period. 

In  addtion  to  ordinary  irregularities  in  the  communication  system  itself,  the 
table  also  reflects  other  prominent  conditions,  such  as: 

CTV  power  failure  (high  voltage  and  ripple) 

Line  switching  between  the  SPS  and  the  Modcomp  processor 
Maintenance  visits 
Tests  from  NDPC. 

In  October  02B  was  affected  by  unwanted  interference  signals  in  the  communi¬ 
cations  system.  Ordinary  operation  resumed  21  October.  After  periods  with 
reduced  reliability,  02C  cable  between  Lillehammer  and  Hornsjtf  was  rerouted 
12  October.  28  October  the  02C  subarray  was  switched  to  the  Modcomp  com¬ 
puter  for  test  recording.  04C  was  affected  week  43,  noise  caused  spikes 
in  the  data  and  B/C  loop  was  found  inoperative.  Attempts  to  replace  the 
modem  and  SLEM  power  supply  did  not  improve  the  operation,  as  errors  were 
still  detected  in  the  modem  test  mode  (C-loop).  A  communications  test  28 
October  proved  satisfactory  operation  without  errors. 

In  November  most  of  the  subarrays  were  affected.  10  November  01 A  and  01 B 
modems  situated  at  the  NDPC  lost  power.  01B  was  also  affected  11-13  November 
due  to  high  ICW  error  rate,  and  25  November  due  to  scheduled  NTA  maintenance 
work.  02B  went  down  6  November  and  remained  down  until  11  November,  probably 
due  to  loss  of  power. 


-  9  - 


After  a  period  in  ordinary  operation  D2C  was  again  switched  to  Modcomp  3  Nover.- 
ber.  18  and  19  November  the  same  suharray  was  affected  by  a  broken  air  cable 
in  the  Mesnalien  area.  Week  48  we  found  the  ICW's  returned  from  040.  The 
reason  was  a  short  circuit  between  two  pairs  of  wires  causing  an  artificial 
loop  to  be  set  up. 

In  December  three  arrays  were  affected.  (''2C,  as  usual  used  in  Modcomp  test 
connnections ,  was  brought  back  in  ordinary  operation  week  49.  030  was  down 
4,5,6  and  7  December  for  unknown  reasons.  A  planned  operation  at  the  NTA 
Elverum  affected  030  and  04C  16  December. 

January.  The  operation  at  Elverum  16  December  where  NTA  phased  in  a  new 
system  had  a  bad  affect  with  respect  to  030,  as  they  were  not  able  to  restore' 
the  original  condition  on  the  communications  system.  From  11  January  02C  was 
used  more  or  less  continuously  in  connection  with  Modcomp  tests.  Errors  were 
observed  in  connection  with  04C  ICW's,  but  the  number  of  errors  was  less  than 
50  and  did  not  affect  or  deteriorate  the  data. 

Also  in  February  02C  was  switched  to  Modcomp  most  of  the  time.  8  February  01B 
was  affected  due  to  low  level  (-40.0  dB).  01 A  was  affected  week  8  due  to  low 
level  toward  the  subarray.  01B  and  03C  were  affected  the  same  week,  causing 
ICW  and  0DW  outage,  respectively. 

March .  As  in  the  previous  months  02C  was  frequently  used  in  connection  with 
Modcomp  tests.  Also  04C  and  03C  were  used,  but  for  shorter  periods.  03C 
was  affected  by  suspected  cable  fault  in  the  Rena  area.  To  restore  the 
B/C  loop  feature,  a  card  was  replaced  in  the  03C  modem. 

Table  II. 4. 2  indicates  distribution  of  outages  with  respect  to  the  individual 
subarrays.  The  table  also  reflects  (by  asterisks)  weeks  when  some  or  all 
subarrays  have  been  affected  simul taneously . 

1 1 . 5  The  ARP  A  Subnetwork 

In  October  and  November  periods  with  reduced  line  quality  were  observed. 
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II .6  The  Terminal  Interface  Message  Processor  (TIP) 

In  connection  with  installation  of  a  power  distribution  panel  the  TIP  lower 
was  removed  for  4-5  hours  10  November.  Week  50  the  NORSAR  conrection  to  the 
ARPA  network  and  TIP  was  terminated.  System  power  was  turned  of  21  December. 

O.A.  Hansen 
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Week/  Subarray/per  cent  outage 


Year 

01A 

01B 

02B 

02  C 

03C 

04  C 

41/81 

36.1 

2.8 

0.7 

0.2 

42 

- 

0.7 

- 

2.1 

- 

1  .2 

43 

0.2 

0.3 

28.7 

- 

0.2 

93.4 

44 

- 

0.4 

- 

28.1 

0.3 

19.2 

45 

- 

- 

42.3 

73.2 

0.1 

- 

46 

4.3 

22.2 

31.2 

100.0 

0.4 

3.1 

47 

- 

- 

- 

100.0 

1.2 

0.4 

48 

0.4 

3.1 

0.4 

100.0 

0.9 

3.3 

49 

- 

2.4 

- 

50.6 

11.9 

0.1 

50 

0.1 

0.1 

0.1 

2.2 

15.6 

0.4 

51 

- 

- 

- 

- 

0.4 

0.4 

52 

1.5 

0.2 

0.4 

0.1 

0.1 

0.3 

53 

- 

- 

- 

- 

- 

2.2 

1/82 

0.1 

- 

- 

2.1 

- 

r  .6 

2 

- 

- 

- 

37.5 

0.6 

«,  .4 

3 

- 

2.2 

- 

95.2 

9.1 

r.5 

4 

0.6 

- 

- 

80.9 

5.5 

0.4 

5 

3.6 

1.8 

- 

100.0 

15.9 

0.3 

6 

46.1 

0.1 

0.1 

8.7 

8.9 

0.7 

7 

- 

0.1 

0.1 

100.0 

7.6 

0.6 

8 

0.6 

2.8 

0.3 

100.0 

57.4 

0.1 

9 

1.0 

1.0 

1.0 

25.0 

71.4 

0.9 

10 

- 

- 

0.1 

- 

71.4 

- 

11 

- 

- 

- 

91.7 

11.4 

1.2 

12 

- 

- 

- 

100.0 

11.1 

0.6 

13 

- 

- 

- 

51  .2 

7.5 

- 

TABLF.  1 1. 4. 2 


-  n  - 

TIT.  IMPROVRMF NTS  AND  MOD  I  FT CAT  TONS 

JII.l  NORSAR  On-line  System  using_  IBM  4331/4341  and  Modcomp  Classic 
The  IBM  VM/SP  operating  system  is  designed  for  fast  terminal  access  to 
many  simultaneous  users,  and  therefore  given  highest  priority  to  handle 
terminal  interrupts.  This  priority  is  huilt  into  the  system  and  may  not 
he  changed.  This  means  that  a  real-time  system,  i.e.,  the  Modcomp  to  IBM 
communication,  running  under  VM/SP  control,  may  he  disturbed  hv  use  of 
a  t  e  rm  i  na  1 . 

For  the  design  of  our  new  DP  system,  this  means  that  we  should  not  allow 
any  terminal  use  on  this  machine,  since  this  may  cause  loss  of  interrupts 
from  the  Modcomp  communication  controller.  For  NORSAR  this  is  no  problem, 
since  all  control  of  the  DP  system  is  done  via  shared  disk  access  between 
the  4331  and  4341  computer.  Moreover,  the  Max4  operating  system  from 
Modcomp  will  later  this  year  support  their  IBM  channel  interface  which 
was  originally  in  our  design  for  the  SPS  substitution.  If  a  change  to  this 
kind  of  interface  is  done,  the  transfer  of  data  may  be  at  a  speed  of  513 
Kbytes  rather  than  230  Kbit  per  second.  The  earlier  reported  loss  of  data  wa 
not  mainly  due  to  the  above  discussed  terminal  access  priority  in  VM/SP,  but 
due  to  the  handling  of  error  checking  in  the  IBM  bisync  protocol  (CROcheck). 
We  can  now  conclude  that  with  correct  CRC  generation  on  Modcomp,  and  no  extra 
users  on  the  IBM  system,  we  have  no  loss  of  data  in  the  transmission  between 
Modcomp  and  IBM. 

We  have  the  following  elements  of  the  on-line  system  in  operation: 
Communication  with  Modcomp.  Receiving  high  rate  data  block,  subarrav  beam  Hat 
and  NORKSS  data  o.k.  Transferring  FS/FRS  data  and  manual  status  to  Modcomp  o. 
No  loss  of  data. 

Recording  of  data  on  disk.  o.k.  We  have  direct  access  vising  a  system  when'  th 
time  of  the  data  determines  which  block  to  use,  i.e.,  given  time  of  data, 
number  of  disks  and  disk  size  we  can  read  the  data  directly  with  no  search. 
The  data  is  available  to  the  4341  by  shared-disk  access.  A  total  of  3b  hours 
may  presently  be  recorded  on  the  allocated  disks.  The  space  allocated  for 
DPX's  may  cover  one  year  of  DPX  data. 
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Detection  processing.  So  far  id  problems  have  been  encountered.  A  few  hours 
of  full  array  detection  test  have  been  performed,  and  a  concluding  test  running 
over  several  days  will  be  performed  in  a  short  time. 

Ta^e  rec_ordjing.  o.k.  6250  BPI  is  used  and  one  tape  will  cover  12  hours  of  NORSAR 
and  NORESS  data.  Presently  a  blocking  of  9  seconds  of  NORSAR  data  is  used. 

The  new  tape  format  will  be  compatible  with  present  high  rate  tapes,  except  for 
a  longtr  blocking  than  the  present  2.0  seconds  block  and  the  presence  of  'NA'~ 
records  (NORESS). 

Calibration  status.  All  software  necessary  to  transmit  any  FS/FRS  is  o.k.  Manual 
status  setting  is  o.k.  All  cal ibrat ion/ status  commands  are  executed  on  the 
4341  computer.  Simulation  of  EOC  alternate  telemetry  commands  is  not  yet 
implemented . 

Experimental  operation  console  (FOC).  '’he  FOC  now  displays  seismometer  data  in 
a  scrolling  manner,  i.e.,  the  traces  w' 11  move  out  of  the  screen  at  one  edge  as 
new  data  is  coming  in  on  th<_  opposite  <  dge  of  the  screen. 

The  hardware  we  have  on  the  4341  system  is  capable  of  displaying  onl  r  ’frozen' 
pictures  (i.e.,  Tektronix  618  screens  attached  to  the  IBM  3277  graphic  attach¬ 
ment).  T  e  plot  software  developed  may  display  any  number  of  short  period  data 
from  NOR:  AR  and/or  NORFSS.  Default  option  will  plot  the  latest  50  seconds 
recorded,  or  any  time  interval  that  is  present  on  the  on-line  disk.  A  program 
for  displaying  long  period  data  is  also  developed.  The  same  program  may  handle 
calibration/status  commands  so  that  you  may  initiate  calibration  of  an  instru¬ 
ment  and  'immediat  •  l y '  see  the  sine  waves  displayed  on  the  graphics  screen. 

All  ROC  commands  end  up  in  FS/FRS  sequences  or  DOF1/DOF2,  and  so  far  we 
have  capabilities  of  initiating  FS/FRS  commands.  The  full  simulation  of 
all  EOC  commands  will  be  implemented  as  needed. 

Note:  Although  we  have  designed  a  system  where  new  terminals  will  be  used 

for  status,  calibration  and  plotting,  we  might  also  consider  attaching 


-li¬ 
the  EOC  to  the  4341  computer  directly.  The  software  Involved  in  using 
the  EOC  as  it  Is  used  in  the  old  system  wi 1  I  in  that  case  need  to  he 
developed,  a  task  which  we  have  not  yet  looked  into. 

So  far  we  may  conclude  that  all  EOC  operations  except  alternate  telemetrv 
commands  have  replacements  in  the  new  system. 

SPS  substitution  (Modcomp).  The  software  developed  on  the  Modcomp  turned  out 
to  have  a  too  high  overhead  to  run  7  subarrays  (7  BSC  lines)  simultaneously. 

The  subcontractor  has  now  reprogrammed  the  ICW  and  ODW-hand  1  i  ng ,  and  7-subarray 
operation  has  been  tested  successfully.  Some  further  refinements  are  being 
worked  on. 

The  new  TOD  unit  had  to  be  returned  to  the  factory  due  to  occasional 
jumps  in  time  (1  second  time  jumps).  Meanwhile  we  are  using  the  oldest 
TOD  for  Modcomp  tests.  This  has  given  some  problems  to  the  communications 
stability,  but  will  be  resolved  as  soon  as  the  new  unit  is  back  from  the 
factory. 

Event  processing.  So  far  we  have  a  system  for  copying  data  intervals  from  on¬ 
line  disk  to  EPX  files  using  the  PPX  list  for  time  intervals.  The  FASTF.P 
program  Is  changed  to  use  disk  files  rather  than  tape  input.  Event  processing 
is  not  fully  automated  but  the  basic  programs  for  event  processing  from  on¬ 
line  data  is  available  now. 

Conclusion.  The  present  software  on  Modcomp  does  perform  the  SPS  tasks,  and 
mav  replace  SPS  as  it  is.  However,  a  few  improvements  should  he  done. 

The  present  software  on  1 HM  4331  does  perform  the  recording,  detection  and 
ca 1 i bra t ion/status  tasks.  Some  minor  changes  are  necessary  to  include 
message,  DPX  and  EPX  transmission  to  Modcomp  and  further  to  AFTAC.  Due  to 
terminal  priority  in  VM/SP  it  is  not  practical  to  use  the  4331  for  other 
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users  together  with  the  on-line  system,  and  a  stan  -alone  program  would 
be  sufficient  to  run  the  new  DP  system.  This  would  make  more  disk  space 
available . 

The  present  software  on  IBM  4341  can  now  perform  the  tape  recording,  on-line 
plotting,  calibratin  and  status  (array  monitoring)  task,  and  event  processing 
may  be  done.  The  latter  needs  some  more  organization  to  use  the  disk  space 
effectively  for  EPX  data. 

The  next  reporting  period  should  give  the  performance  of  and  experiences 
with  the  new  system. 


J .  Fyen 


III .2  Field  Instrumentation  and  facilities 

Refer  to  Table  III. 2.1  for  changes  of  the  array  stations  up  to  the  end  of 
May  1982.  An  analog  line  from  06C  CTV  was  operational  as  of  8  October  1981 
carrying  eight  of  the  NORESS  channels  (see  Table  III. 4.1  in  Torstveit,  1981). 
The  channels  are  digitized  and  recorded  by  the  MODCOMP  computer.  The  Landrover 
was  replaced  with  a  new  Mercedes  300D  at  the  end  of  the  reporting  period  and 
fully  operational  in  May. 


A.  Kr.  Nilsen 


Reference 
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Status  of  NORSAR  SP  instruments  recorded  on  data  tape. 
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IV.  FIELD  MAINTENANCE  ACT  TV  TTY 

In  this  section  a  very  brief  review  of  the  field  maintenance  is  given.  Ther.- 
'  .  little  change  from  previous  periods  and  the  error  rate  and  maintenance 
activity  are  normal.  Two  new  stations  are  operative  in  the  subarray  12B 
area  and  radio  transmitted  to  NOREAR  field  maintenance  center  at  Stange 
and  recorded  on  a  Kinemetrics  PDR-2  digital  seismic  recorder. 

Maintenance  Visits 

Table  IV . L  gives  the  number  of  visits  to  the  NORSAR  subarrays  during  the 
reporting  perod.  The  average  number  of  visits  to  each  subarray  is  2.4,  not 
including  the  work  with  the  new  stations  in  the  02B  area. 


Subarrays  01A  01B  02B  02C  0_3C  04C  __  06C/NQRESS  Total 

No .  of  Visits  __  5^ _ j4__ _ _ _ 1_ _ 2_ _ 1 _ _ 4_  _ _ _ 11  _ 

TABLE  IV. 1 

Number  of  visits  to  the  NORSAR  subarrays  including  NORESS 
in  the  period  1  Oct  1981  -  31  Mar  1982 
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Preventi'-e  Maintenance  Projects 

The  preventive  maintenance  work  of  the  NORSAR  array  instrumentation  and  facil¬ 
ities  listed  in  Table  TV'. 2.  The  adjustments  are  corrections  of  characteristics 
within  the  tolerance  limits. 


Unit 


Action 


No .  of 
Actions 


Seismometer 


MP  adjust  (in  field) 

pp 


14 

1 


Line  Termination  Adjustment  of  channel  gain  (SP) 
Amplifier  DC  offset  (SP) 


Emergency  Power  Battery  and  charger  check 


Cleaning  of  CTV 


1 


TABLE  IV. 2 


Preventive  maintenance  work  in  the  period 
1  October  1981  -  31  March  1982 
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Corrective  Maintenance 

The  corrective  maintenance  comprising  the  required  adjustment  and  replacements 
is  given  in  Table  TV'. 2. 


Unit 

Characteristics 

_ S P  _ 

Repl .  Ad j. 

_ LP 

Repl .  Ad  j . 

Seismometer 

FP  (in  field) 

14 

MP  (in  field) 

6 

MP/FP  (at  NDPC) 

43 

Ren 

o 

Seismometer  Ampl. 

Gain 

1 

RA/5,  Ithaco 

Protection  card 

1 

Taperpin  block 

1 

Line  Termination 

Gain 

Amplifier 

DC  offset 

1 

EPU  1 


SLEM 

Mux  Spikes  1 


TABLE  IV. 3 


Total  number  of  required  adjustments  and  replacements  of  NORSAR  field 
equipment  in  the  period  1  October  1981  -  31  March  1982. 
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Powe r  Breaks ,  Cable  Breakages ,  Commun i c at i on  s  Fault s 

One  power  break  required  action  of  our  field  technicians,  and  brol en  cables 
to  01B  04,05  were  repaired.  Two  cominunicat ion  faults  were  correcte  1  by  our 
field  crew. 


Array  Status 

As  of  31  March  1982  the  following  channels  had  out-of-tolerance  conditions: 
01A03 ;  02B04;  03C06,V;  04C  V;  NORF.SS  06,09;  06C  V,NS. 

Channels  with  nonstandard  conditions  are: 

01A  04 

01A  06 

02  B  06 
04C  06 

06C  01,03-06 


Attenuated  30  dB 
8  Hz  Filter 

__M  _ 

NORESS  (02  is  placed  in  LPV) 


A.Kr.  Nilsen 
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ABBREVIATIONS 


CTV 

- 

Central  Terminal  Vault 

DC 

- 

Direct  current 

EPU 

- 

External  Power  Unit 

FP 

- 

Free  period 

LP 

- 

Tong  period 

MP 

- 

Mass  position 

MUX 

- 

Multiplexer 

NDPC 

- 

NORSAR  Data  Processing  Center 

NORESS 

- 

NORSAR  Experimental  Small-Aperture  Subarray 

RCD 

- 

Remote  centering  device 

SLEM 

- 

Seismic  short  and  long  period  electronics  module 

SA 

- 

Subarray 

SP 

- 

Short  period 

24 
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VI.  SUMMARY  OK  TKCHMCAI.  RK PORTS/ PARKRS  PR F. PARK!) 

VI. 1  Seismic  source  spectra  and  moment  tensors 

The  quantification  of  seismic  sources  ( e .  g .  ,  scalar  moment  or  magnitude), 
the  determination  of  source  size  (e.g.,  corner  frequency  methods),  and 
the  determination  of  radiation  pattern  (e.g.,  fault  plane  solutions), 
are  procedures  which  usuallv  proceed  separately  and  use  different  parts 
of  the  seismic  spectrum.  Of  course,  the  various  source  parameters  are 
not  independent  of  one  another,  and  they  are  often  related  through  semi- 
empirical  rules  (Kanamori  and  Anderson,  1975).  There  would  be  obvious 
advantages  if  source  parameters  could  be  determined  s imul taneousl v.  The 
seismic  moment  tensor  represents  both  scalar  moment  and  radiation  pattern, 
but  applications  are  restricted  to  what  may  be  regarded  as  point  sources. 
In  principle  it  is  possible  to  remove  this  restriction  by  extending  tbe 
representation  to  moment  tensors  of  higher  degree  (Backus,  1977),  and 
as  described  in  the  previous  semiannual  report,  in  practice  it  is 
possible  to  estimate  moment  tensors  up  to  degree  two  (Doornbos,  19R2a). 

At  this  stage,  however,  a  long-period  approximation  is  still  implied, 
and  spectral  information  above  the  corner  frequency  cannot  be  used  in 
the  present  procedure. 

In  a  more  recent  attempt  to  reconcile  the  information  from  long-  and 
short-period  data  (Doornbos,  1982b)  it  was  shown  that  when  the  source 
spectrum  is  expanded  in  powers  of  frequency,  the  coefficient  of  the 
n'th  power  is  a  linear  function  of  the  n'th  moment  tensor: 

F(£,w)  =  1-iio£TF(  |  )-it°2£T^(2)5.+*  *  *  ( 1  ) 

where  F(i)  and  _F(2)  are  the  moment  tensors  of  degree  1  and  2,  respectively 
and  C  is  a  4-dimensional  slowness  vector: 

£  =  (Y,/c,  Y2/c,  Yq/c,l)T 

where  c  is  wave  velocity  and  a  direction  cosine  of  the  wave.  It  should 
he  noted  that  equation  (1)  gives  the  moment  tensor  expansion  of  a  scalar 
source  function  f(£,l)  which  is  relative  to  the  moment  tensor  density  of 
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Ajk(i,T>  =  Mjkf(S,T>  (2) 

and  Mj^  Is  the  moment  tensor  of  degree  zero. 

The  form  (1)  explicitly  shows  the  low-frequency  approximat ' on  by  moments 
of  low  degree,  and  suggests  suitable  ways  of  extrapolation  In  e <trapolating 
to  higher  frequencies  the  implied  spectral  assumption  rest  -icts  the  source 
model  to  a  particular  class,  and  in  Fig.  VI. 1.1  several  possibilities  are 
compared  on  the  basis  of  a  typical  triangular  earthquake  source  pulse.  The 
models  included  are  among  others  oi-square: 

HU,u.)  =  |l+iw2£TJ(2)_£}_1  exp{-im?'IV(1)}  (1) 

and  Gaussian: 

N(c,w)  =  exp{“iu2c’^F(2)£  ~  lwiTi'(i)}  (M 

Here,  2)  is  measured  with  respect  to  the  source's  'center  of  gravity' 
in  space  and  time  (Cq.t'q),  for  which  F(^)  =  0.  The  model  has  P0  para¬ 
meters,  involving  moments  of  degree  zero,  one  and  two.  Model-dependent 
constraints  can  reduce  this  number  to  that  employed  in  other  methods. 

The  models  lead  to  simple  expressions  for  the  displacement  field  and  the 
total  radiated  seismic  energy,  and  they  can  be  made  to  satisfy  certain 
general  properties  of  observed  far-field  spectra,  including  the  corner 
frequency  shift  of  P  waves  with  respect  to  S  waves. 

Moments  of  degree  zero  represent  the  final  static  source  parameters, 
which  are  usually  obtained  in  a  point  source  approximation  from  long- 
period  data.  A  mislocation  of  the  source  gives  rise  to  moments  of  degree 
one.  They  determine  the  phase  spectrum  in  the  spectral  model,  and  their 
determination  is  equivalent  to  the  class’cal  source  location  problem  by 
travel  time  analysis.  The  source's  spatial  and  temporal  extent  give 
rise  to  moments  of  degree  two.  They  control  the  spectral  bandwidth  and 
its  variation  with  take-off  angle  and  wave  velocity,  and  their  determination 
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corner  frequency  methods  of  determining  source  dimensions.  Although 
solutions  for  the  moments  of  different  degree  are  not  independent,  it 
should  often  be  possible  to  obtain  initial  estimates  of  the  moments  of 
degree  zero  (linear  inversion),  one  and  two  (linear  or  nonlinear  inversion 
consecutively,  from  long-  and  shortperiod  data  or  more  precisely,  from 
data  in  at  least  two  frequency  bands. 

Inversion  methods  based  on  this  model  are  being  developed  to  obtain 
estimates  of  scalar  moment,  seismic  energy,  source  size  and  stress  drop, 
from  long-  and  short-period  SRO/ASRO  data.  It  is  also  found  that  the 
following  points  need  further  study:  (1)  The  presently  employed  digital 
SRO  response  can  be  small  in  crucial  parts  of  the  excitation  spectrum  of 
commonly  observed  sources  (m^,  roughly  in  the  range  5-6).  (2)  Frequency 
dependence  of  Q  appears  to  be  needed  in  source  retrieval  from  long- 
and  short-period  data.  (3)  Short-period  amplitude  anomalies  can  signifi¬ 
cantly  affect  the  source  estimates.  Points  (2)  and  (3)  corroborate  the 
results  of  others. 

D.J.  Doornbos 
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Fig.  VI. 1.1  Spectral  approximations  to  triangular  source  pulse  (see  inset) 
with  amplitude  spectrum  F(w)  and  corner  frequency  u  . 

H(w):  w-square  model  (eq.  O));  N(u>):  Gaussian  mode'l  (eq.  ( /, ) ) 
P(w)  and  K(w)  are  a  parabolic  and  a  rational  approximation, 
not  mentioned  in  the  text. 
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VI. 2  A  processing  package  for  on-line  analysis  of  data  from  small- 
aperture  arrays 

A  Regional  'ON-line'  Array  Processing  Package  (RONAPP)  has  been  developed 
and  tested  on  off-line  data  from  the  smal l-aperture  NOKKSS  array  at 
■VORSAR.  This  array  consists  of  12  sensors  spaced  with  distances  up  to 
about  2  km,  and  the  processing  package  automatically  detects  and  locates 
regional  seismic  events  in  the  distance  range  from  a  few  kilometers  to 
the  maximum  propagation  range  of  regional  S-tvpe  waves  (Sn  or  I,g).  The 
package  consists  of  a  conventional  STA/LTA  detector,  a  phase  identifica¬ 
tion  procedure  based  on  phase  velocity,  and  a  location  procedure  based 
on  observed  travel  time  differences  and  a  common  azimuth  between  the 
observed  primary  and  secondary  phases.  In  the  following  we  give 
some  details  on  the  detection  and  detection  processing  procedures  and  a 
description  is  given  of  the  phase  association  scheme  adopted.  Finally 
some  processing  results  are  shown  and  discussed. 

Detection  and  dot < “ c tio n  pro e e  s_s_i_n g 

In  the  present  version,  the  detector  works  on  either  a  vertical  beam  or 
a  single  channel  and  with  a  number  of  filters.  For  each  filter,  STA/I.TA 
is  computed  and  compared  successively  against  two  different  thresholds. 
Three  blocks  of  data  (presently  with  A  sec  in  each)  are  always  kept  in 
memory,  with  the  detection  performed  in  the  middle  block.  If  no  detection 
is  declared,  the  first  block  is  dropped,  a  npw  one  is  brought  in,  and  the 
procedure  is  repeated.  In  case  of  one  or  more  detections  in  the  current 
block,  data  within  a  certain  time  window  around  each  detection  are  written 
on  disk,  estimates  of  refined  arrival  time  and  signal  frequency  are  ob¬ 
tained,  and  a  wavpnumber  analysis  is  performed  following  a  proper  pre- 
filtering.  This  gives  estimates  of  phase  velocity,  azimuth  and  power  ror 
each  detection,  which  together  with  arrival  time  and  frequency  constitutes 
the  essential  part  of  a  detection  record.  The  detection  procedure  is 
continuous  and  not  dependent  on  the  block  structure.  A  block  diagram 
outlining  the  detection  and  detection  processing  procedures  is  shown  in 
Fig.  VI. 2.1. 
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Phase  association  and  event  location 

Following  each  processed  detection,  an  event  location  is  performed  if 
the  last  detection  is  identified  as  an  S-type  wave  (phase  velocity  below 
6  km/s),  and  also  a  P-wave  (phase  velocity  above  6  km/s)  has  been  detected 
within  an  appropriate  azimuth  range.  Distance  is  taken  from  an  S-P  travel 
time  table.  An  intricate  procedure  is  developed  to  allow  for  several 
detections  within  the  P  and  S  coda,  to  keep  track  of  'used'  and  'unused' 
detections,  and  to  identify  locations  that  override  previous  ones. 

The  identification  of  phase  type  according  to  phase  velocity  being  above 
or  below  6  km/sec  is  justified  from  'manual'  wavenumber  analysis  of  NORESS 
data  from  well-identified  phases  (Mykkeltveit  &  Ringdal,  1981). 

The  automatic  phase  association  and  location  procedure  assumes  in  its 
present  version  that  the  largest  secondary  arrival  detected  is  an  Lg 
phase.  The  block  diagram  in  Fig.  VI. 2. 2  gives  the  details. 

The  main  problem  in  the  location  procedure,  however,  is  the  correct 
identification  of  the  type  of  secondary  arrival  (in  practice  Sn  or  Lg) 
which  is  associated  with  P.  It  is  not  possible  from  our  wavenumber  analysis 
to  distinguish  between  Sn  and  Lg  on  the  basis  of  phase  velocities  alone. 
Unless  other  objective  criteria  can  be  found,  it  will  be  necessary  to 
invoke  region-dependent  information  on  the  general  occurrence  of  secondary 
seismic  phases.  Regional  events  within  12-15°  recorded  by  NORESS  are  domi¬ 
nated  by  strong  Lg  waves,  propagating  at  a  fairly  constant  group  velocity 
of  3.5  km/sec.  Exceptions  are  events  in  Central  Europe,  England,  Scotland 
and  western  part  of  the  North  Sea,  from  which  Sn  is  the  dominating  phase 
and  Lg  is  very  weak  if  at  all  seen.  The  strong  attenuation  of  Lg  is 
probably  an  effect  of  larterally  varying  structures  along  paths  to  NORESS 
from  these  regions.  Events  from  or  close  to  the  North  Atlantic  Ridge 
show  the  same  characteristics.  A  detailed  Investigation  of  regional 
propagation  characteristics  within  15°  of  ,«0RSAR  and  correlation  with 
gross  geological  features  is  In  progress. 
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Experience  with  the  processing  package 

The  processing  package  has  been  tested  on  a  number  of  data  intervals  from  th 
12-sensor  NORESS  array.  Our  experience  so  far  can  be  summarised  as  follows: 
The  on-line  automatic  processing  package  provides  location  results  as  good 
as  those  previously  reported  by  Mvkkeltveit  and  Ringdal  (1981)  for  manual 
analysis  of  NORESS  data,  provided: 

a)  proper  positioning  of  time  window  for  f-k  analysis 

b)  the  automatically  determined  frequency  for  the  f-k  analysis  must  corre¬ 
spond  to  peak  in  signal  spectrum 

c)  data  to  be  subjected  to  f-k  analysis  must  be  properly  prefiltered 

d)  assumptions  inherent  in  the  automatic  phase  association  are  '"lid. 

Points  b)  and  c)  ahove  ensure  reliable  phase  velocity  and  azimuth  results 
from  the  f-k  analysis,  while  d)  is  an  obvious  requirement  for  proper 
location.  Fig.  VI. 2. 3  serves  as  an  illustration  of  the  importance  of  the 
positioning  of  the  time  window  for  f-k  analysis.  This  figure  gives  the 
results  from  a  sliding  window  analysis,  where  phase  velocity  and  azimuth 
results  are  plot'ed  for  windows  of  length  2.3  sec  shifted  0.6  sec  at  a  time. 
Each  time  w'ndow  has  been  subjected  to  f-k  analysis  for  9  frequencies  in 
the  range  3-5  Hz  and  results  plotted  (at  the  time  interval  midpoint)  cor¬ 
respond  to  the  frequency  with  maximum  power. 

Fig.  VI. 2. 3  shows  that  azimuth  and  phase  velocity  estimates  for  this  high 
signal-to-noise-rat io  event  are  comparatively  stable  around  the  onset  times 
for  the  main  phases.  Azimuth  values  derived  from  time  intervals  in  the  codas 
on  the  other  hand,  fluctuate  around  the  expected  value.  For  this  event, 
however,  all  P  and  Eg  coda  detections  (two  for  each  phase)  are  associated 
with  lower  power  than  the  phase  onset  detections,  which  are  marked  in  tin1 
figure.  Consequen t 1 y ,  our  association  and  location  procedures  work  well 
in  this  case. 


Epicentral  distance  is  determined  from  the  first  P  and  the  strongest  hg- 
type  detection  via  the  time  difference  between  these  two  detections.  A  good 
estimate  of  distance  relies  upon  the  assumption  that  the  amplitude  maximum 
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in  the  Lg  wave  train  can  be  associated  with  a  constant  group  velocity 
of  3.5  km/s.  This  assumption  has  been  justified  by  the  testing  performed 
so  far. 

■'ig.  VI. 2. 4  shows  the  data  for  the  12  NORF.SS  channels  and  the  vertical 
oeam  for  a  presumed  mining  explosion  located  by  the  Finnish  network  to 
the  Fi  land-U . S . S . R.  border  at  a  latitutde  of  about  61°N.  The  arrows 
indica  e  the  four  detections  by  our  processing  package  and  results  from 
f-k  analysis.  Besides  the  P  and  Sn  detections,  there  are  two  Lg  detections; 
the  f-k  plot  shown  corresponds  to  the  strongest  one.  This  example  shows 
that  our  detection  and  location  procedure  works  also  for  a  fairly  modest 
signal-to-noise-ratio  event. 

Although  more  refinement  and  testing  are  still  needed,  we  hive  finalized 
the  basic  design  for  a  package  that,  when  finished,  should  >e  suitable  for 
implementation  within  the  framework  of  a  low-power,  microconputer-based 
system. 


S.  Mykkeltveit 
H.  Bungum 
F.  Ringdal 
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Y_I-3_  Seismic  array  configuration  optimization 

The  seismological  concept  of  arrays  emerged  from  the  nuclear  test  ban 
negotiations  in  Geneva  in  1958  when  the  need  for  improved  capability 
to  study  weak  seismic  events  was  clearly  recognized.  Over  the  years 
many  kinds  of  array  systems  have  been  installed  in  many  countries,  but 
LASA  and  NORSAR  remained  the  unique  ones;  unique  in  terms  of  their  sizes 
(apertures  of  the  order  of  200  and  100  km,  and  number  of  instruments 
once  588  and  198,  respectively),  and  because  dedicated  data  centers 
were  made  an  integral  part  of  these  arrays. 

In  Norway  a  miniarray  termed  NORESS  (NORSAR  Experimental  Sma1 1  Subarray) 
has  recently  been  installed  within  the  NORSAR  array  itself  and  is  shown  in 
Fig.  VI. 3.1.  The  motivation  for  this  undertaking  is  to  provide  data  for 
optimum  design  of  a  complete  array  system  for  surveillance  of  events 
at  regional  distances.  This  is  the  research  topic  dealt  with  in  this 
section. 

SignaJL  and  noise  structure 

An  essential  aspect  of  the  array  design  is  a  priori  assumptions  regarding 
the  statistical  structure  of  signal  and  noise.  Generally  speaking,  the 
array  response  when  'steered  at'  the  signal,  should  be  small  at  wavenumbers 
characteristic  of  the  noise.  Intersensor  spacings  would  then  correspond 
to  maximum  correlation  for  the  signal  and  a  minimum  for  the  noise. 

A  common  assumption  has  been  that  of  equal  signals  and  uncorrelated 
noise  beyond  a  critical  space  lag,  and  for  sufficiently  large  intersensor 
spacing  the  array  design  criteria  are  trivially  met.  However,  observations 
usually  contradict  these  assumptions,  and  our  optimization  procedure 
will  be  based  on  actually  observed  signal  and  noise  correlations  using 
NORESS  recordings. 

Signal  correlations 

A  major  objective  of  our  work  was  to  derive  an  optimum  layout  for  detect¬ 
ing  events  at  regional  distances,  so  signal  correlations  w^re  measured  for 
phases  identified  as  Pg ,  Pn  and  Lg  for  the  5  events  listed  in  Table  VI. 1.1, 
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all  within  12°  of  NOiESS.  A  record  section  of  Event  3  is  displayed  in 
Fig.  VI. 3. 2,  signal  correlations  in  Fig.  VI. 3. 3.  We  see  that  the  Pg 
phas »  exhibits  high  correlation  values  for  all  frequency  bands  and 
sens >r  separations  displayed.  The  Lg  phase  on  the  other  hand  correlates 
rather  poorly  in  the  higher  frequency  ranges.  The  results  for  the  Pn 
phase  show  increasing  correlation  with  frequency,  which  reflects  a 
corresponding  peak  in  the  signal  spectrum  around  4  Hz. 

Noise  correlation 

Selected  noise  records,  including  both  summer  and  winter  conditions, 
were  subjected  to  correlation  analysis,  and  the  results  obtained  are 
also  given  in  Fig.  VI. 3. 3.  The  most  notable  feature  here  is  that  the 
correlation  curves  consistently  have  a  negative  minimum  before  tending 
to  zero,  which  is  rated  significant.  This  observational  fact,  easily 
missed  when  noi:  e  similarity  is  expressed  via  the  coherency  measure, 
is  consistent  with  a  model  of  propagating  noise,  as  demonstrated 
by  Mykkeltveit  et  al  (1982). 

In  Fig.  VI. 3. 4  we  show  the  theoretical  correlation  function  for  noise 
with  a  uniform  wavenumber  spectrum  in  the  range  (2tt)_1  to  5»(2tt)~J  c/km. 
This  wavenumber  range  would  be  appropriate  for  typical  Rayleigh  wave 
phase  velocities  (3-4  km/s)  in  a  frequency  range  roughly  0.6  to  3  Hz, 
or  typical  body  wave  phase  velocities  (above  6  km/s)  at  frequency 
above  1  Hz.  The  curve  is  reasonably  close  to  the  observed  curve  which  is 
taken  from  Fig.  VI. 3. 3  and  testifies  to  the  importance  of  propagating 
noise  at  the  Rayleigh  wave  velocities,  even  at  relatively  high  frequencies. 
This  is  in  accord  with  the  results  of  wavenumber  analysis,  although  higher 
phase  velocities  typical  of  body  waves  are  also  observed. 

The  optimization  procedure 

Let  signal  to  noise  ratio  gain  by  beamforming  be  expressed  by 


C2  =  l  WJW.CJ,/},  WjWjPjj 

ij  ij 


(1) 


where  c-jj  Is  the  signal  correlation  between  sensors  i  and  i, 
similarly  describes  the  noise  correlation  and  wj  are  sensor  weights. 

In  the  common  beamforming  practice  at  arrays,  w^=l.  Then,  since  cj^ 

=  =  1,  equation  (1)  can  he  rewritten  in  the  more  common  form: 

l+(N-l)c 

C,2  = _  (2) 

l+(N-l)p 

where  c,  p  describe  average  signal  and  noise  correlations  for  an  N  sensor 
array.  Alternatively,  it  is  possible  to  use  weighted  beamforming  (e.g., 
Christof fersson  and  Husebye,  1974),  and  the  optimum  gain  function  for  noise 
suppression  is  simply  the  reciprocal  maximum  likelihood  estimate  of  noise 
power  (c.f.  Boverie  and  Gregg,  1971): 

C’2  =  l  Pi] 
ij  ' 

where  form  the  inverse  correlation  matrix  of  the  noise.  Our  reasons 

for  optimizing  G  rather  than  G'  are  that  (i)  G  with  unit  weights  reflects 
the  common  beamforming  practice  at  arrays,  (ii)  G  contains  both  noise  and 
signal,  (iii)  G  is  more  stable  with  respect  to  variations  in  the  estimate 
of  noise  correlation. 

The  gain  is  a  nonlinear  function  of  the  sensor  coordinates,  and  for 
more  than  three  sensors  we  obtained  a  maximum  of  G  by  applying  a  standard 
rapid  descent  method  (Fletcher  and  Powell,  1963)  to  analytic  approximations 
of  the  correlation  functions  in  Fig.  VI. 3. 3.  As  might  be  expected,  the 
optimum  configuration  is  strongly  determined  by  the  position  of  the 
noise  correlation  minimum,  but  for  relatively  poorly  correlating  signals 
like  Lg,  the  slope  of  the  signal  correlation  curves  is  also  important. 

Since  these  features  are  frequency  dependent  (see  Fig.  VI. 3. 3),  different 
optimum  geometries  will  result  for  different  frequency  bands.  A  compromise 
is  to  optimize  a  combination  of  the  gains  in  different  frequency  bands. 
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Figs.  VI. 3. 5  and  VI. 3. 6  show  examples  of  the  gain  averaged  over  the 
5  frequency  bands  of  Fig.  VI. 3. 3,  where  the  optimization  was  done  se¬ 
quentially  from  3  to  20  sensors.  Whereas  the  gain  for  Lg  is  near  the 
standard  factor  /N  (in  amplitude),  which  results  from  the  assumption 
of  identical  signals  and  uncorrelated  noise,  the  gain  for  Pg  is  system¬ 
atically  higher. 

Our  final  configuration  with  20  sensors  for  Pg  gives  a  gain  4  dB  above  /N. 

As  demonstrated  in  Fig.  VI. 3. 7,  together  with  the  gains  in  individual 
frequency  bands.  For  example,  in  Fig.  VI. 3. 7  we  have  included  the  optimum 
geometry  for  the  frequency  band  1. 6-4.0  Hz,  with  clearly  reduced  intersensor 
spacings  as  compared  to  the  averaged  geometry  for  the  5  frequency  bands. 
Corresponding  results  for  Lg  are  given  in  Fig.  VI. 3. 8;  here  the  optimum 
geometry  for  frequency  band  1. 6-4.0  Hz  has  larger  dimension  tl  an  the 
averaged  geometry. 

Obviously,  the  optimum  3  sensor  array  should  form  an  equilateral  triangle, 
with  intersensor  spacing  close  to  the  noise  correlation  minimum.  From 
Figs.  VI. 3. 5  and  VI. 3. 7  it  is  clear  that  for  Pg  the  same  optimum  intersensor 
spacing  dominates  the  array  geometry  for  a  larger  number  of  sensors.  An 
illustrative  example  is  the  optimum  geometry  for  20  sensors  (Fig. 

11.3.7),  with  two  pairs  of  sensors  near  the  center  of  the  array.  This 
is  interpreted  as  indicating  the  importance  of  these  positions,  with  a 
maximum  number  of  other  sensors  at  optimum  distance.  For  Lg  (Figs.  VI. 3. 6 
and  VI. 3. 8),  similar  observations  can  be  made  with  regard  to  number  of 
sensors  and  array  dimension,  the  principal  difference  being  that,  due 
to  reduced  signal  correlations,  intersensor  spacings  and  array  dimension 
are  also  reduced.  More  generally,  the  trade-off  be .ween  the  effect  of 
increased  number  of  sensors  and  increased  array  dimension  should  be 
considered  for  limiting  the  final  number  of  sensors  in  the  regional 
array.  Another  important  point  is  that  the  array  performance  inevitably 
would  depend  on  phase  type  and  frequency,  but  'worst  case'  events  .an 
be  considerably  improved  by  flexible  weighting  schemes. 


Although  the  main  scope  of  our  study  has  been  array  design  for  sur¬ 
veillance  of  regional  events,  the  methodology  applied  is  not  limited 
to  this  task.  For  example,  our  array  concept  is  also  the  likely  basic 
subarray  unit  in  a  large  aperture  ( ' teleseismic ’ )  array,  in  which  case 
the  'optimum'  main/side  lobe  strategy  for  deployment  of  subarray  units 
as  once  used  for  LASA  and  NORSAR  would  be  adequate.  A  more  comprehensive 
presentation  of  the  work  synthesized  above  is  given  by  Mykkeltveit  et  al 
(1982). 

S.  Mykkeltveit 

K.  Astebtfl 

D. J.  Doornbos 

E. S.  Husebye 
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Date 

Origin  time 

Location 

Local  Magnitude 

06 

Nov 

1980 

14 .53.02 

59 . 5°N 

10.7°E 

2.1 

25 

Nov 

1980 

02.39.49 

58.4°N 

13.7°E 

2.4 

29 

Nov 

1980 

^  '2.16 

5 1 . 2°N 

18.5°E 

3.5 

26 

Feb 

1981 

17.43.53 

60.3°N 

15.9  °E 

2.1 

01 

Mar 

1981 

05.08.16 

62.8°N 

6.2°E 

2.7 

Table  VI. 3.1 

Local  events  used  in  this  study 
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Fig.  VI. 3. 2  NORF.SS  records  (channels  1-6  in  Fig.  VI. 3.1)  for  event  5  in 
Table  VI. 3.1.  For  each  phase,  correlation  measurements  are 
made  from  2  sec  of  data  around  the  phase  maximum.  All  traces 
are  normalized,  numbers  to  the  right  give  scaling  factors.  The 
gap  between  Pg  and  Lg  is  29  sec. 
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Fig.  VI. 3. 3  Correlation  vs  distance  for  Pg,  Pn,  Lg  and  noise  at  NORESS 

for  the  five  frequency  ranges  indicated.  Each  curve  is  based 
on  measurements  from  66  combinations  of  sensor  pairs. 


..  Theoretical  model 
—  —  Observed  noise  coneiafoon 


0^  2000 
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Fig.  VI. 3. 4  Theoretical  correlation  function  for  noise  with  a  uniform 

wavenumber  spectrum  in  the  range  (2tt)-^  to  5 ( 2tt ) ~ ^  c/km.  The 
observed  noise  correlation  function  (dashed  line)  is  reproduced 
from  Fig.  Vi. 3. 3  (frequency  range  1. 6-4.0  Hz). 
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Fig.  VI. 3. 5  Gain  for  the  Pg  phase  as  a  function  of  number  of  sensors  for 

geometries  derived  from  our  optimization  procedure.  The  starting 
geometry  for  3  sensors  along  with  resulting  geometries  for  11 
and  16  sensors  are  shown.  The  gain  is  an  average  over  the  5 
frequency  bands  in  Fig.  VI. 3. 3. 


5  10  15  20 

NUMBER  OF  SENSORS 


Fig.  VT.3.6  Same  as  Fig.  VT.3.S,  but  for  the  Lg  phase. 
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Pg  gains  for  an  optimized  20-eiement  array.  The  lower  right 
geometry  is  optimum  when  the  gain  function  is  expressed  as  the 
sum  of  5  individual  gain  functions  for  the  frequency  bands 
shown  in  Fig.  VI. 3. 3.  This  geometry  gives  different  gain  for 
each  band,  as  indicated  by  the  'single  frequency  gains’.  The 
geometry  derived  from  a  'mulitple'  gain  function  is  different 
from  that  based  on  a  single  frequency  band.  This  is  demonstrated 
by  reoptimizing  the  lower  right  configuration  for  the  frequency 
band  centered  at  2.8  Hz  resulting  in  a  net  gain  of  4  dB  and  a 
new  geometry  as  shown.  The  exceptionally  large  gain  for  the 
frequency  band  centered  at  2.2  Hz  reflects  the  pronounced 
negative  noise  correlation  observed  here  (see  Fig.  VI. 3. 3). 
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Fig.  VI. 3. 8  Lg  gains  for  an  optimized  20-element  array.  The  geometry  to 
the  right  is  optimum  when  the  gain  function  is  expressed  as 
the  sum  of  3  individual  gain  functions  in  the  frequency  bands 
given  in  Fig.  VI. 3. 3.  This  geometry  gives  different  gains  for 
each  band,  as  Indicated  by  the  'single  frequency  gains'.  The 
optimum  geometry  derived  from  a  'multiple'  gain  function  is 
different  from  that  based  on  a  single  frequency  band.  This  is 
demonstrated  by  reoptimizing  the  lower  right  configuration 
for  the  frequency  band  centered  at  2.8  resulting  in  a  net 
gain  of  nearly  8  dB  and  a  new  geometry  as  shown.  The  higher 
frequency  ranges  are  seen  to  dominate  the  'average'  geometry. 


48 


VI. 4  Upper  mant_le  seismic  heterogeneities  beneath  Fennoscand  i_a 
Three-dimensional  seismic  mapping  of  the  upper  mantle  hen >ath  Fennoscandia 
(Baltic  Shield)  using  an  ACH-type  of  inversion  technique  in  comb  nation 
with  P-vave  travel  time  residual  observations  from  the  local  seismograph 
network  gave  the  following  results.  The  central  parts  of  the  Baltic  Shield 
are  characterized  by  relatively  high  seismic  velocities  down  to  approximately 
300  km.  Those  parts  of  the  shield  most  affected  by  the  Caledonide  orogeny 
exhibit  relatively  low  velocities  particularly  in  the  uppermost  100  km 
depth  interval.  The  lower  part  of  the  upper  mantle  (300-600  km)  does  not 
exhibit  pronounced  seismic  velocity  anomalies  and  in  this  respect  is  in 
contrast  to  results  from  similar  studies  in  regions  subjected  to  neotectonic 
processes  like  parts  of  central  and  southeastern  Europe.  The  seismic 
anomaly  pattern  in  the  presumed  thickened  lithosphere  is  in  quantitative 
agreement  with  similar  ones  derived  from  surface  wave  dispersion  analysis 
and  inversion  of  electrical  measurements.  The  general  orientation  of 
these  anomalies  coincides  with  that  of  the  glacial  uplift.  \  comprehensive 
description  of  this  work  can  be  found  in  a  recent  paper  by  Husebye  and 
Hovland  (1982). 


E.S.  Husebye 

J.  Hovland,  Univ.  of  Oslo 
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Husebye,  E.S.  and  J.  Hovland  (1982):  On  upper  mantle  seismic  heterogeneities 
beneath  Fennoscandia.  Tectonophysics ,  in  press. 


VI.5_  Upper  mantle  heterogeneities  beneath  Eastern  Europe 
P-wave  travel  time  residuals  for  seismograph  stations  in  Eastern  Europe 
as  reported  by  ISC  for  the  years  1964-77  were  used  for  constructing  a 
seismic  image  of  upper  mantle  heterogeneities  in  the  network  region. 

For  the  depth  range  0-100  km,  dominant  tectonic  features  like  the  Pannonian 
Basin  and  the  Aegean  Sea  and  Western  Turkey  correlate  well  with  pronounced 
velocity  lows  which  appear  to  extend  down  to  a  TOO  km  depth.  The  velocity 
anomaly  patterns  in  the  depth  intervals  300-500  km  and  500-600  km  are 
broadly  similar  but  quite  different  from  those  of  shallower  depths.  The 
observed  seismic  heterogeneities  are  tniefly  discussed  in  terms  of  large- 
scale  tectonic  and  geophysical  (heat-flow)  characteristics  of  Eastern 
Europe.  A  comprehensive  description  of  this  work  can  be  found  in  a  recent 
paper  by  Hovland  and  Husebye  (1982). 


J.  Hovland,  Univ.  of  Oslo 
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Hovland,  J.  and  E.S.  Husebye  (1982):  Upper  mantle  heterogeneities  beneath 
Eastern  Europe.  Tectonophys ics ,  in  press. 


VI. 6 A  North  Sea  -  Southern  Norway  seismic  crustal  profile 

Introduction 

During  the  last  decade  numerous  seismic  profiling  investigations  have  been 
carried  out  in  Fennoscandia  and  a  relatively  large  number  of  crustal 
models  have  been  derived.  There  is  considerable  inconsistency  concerning 
the  number  of  crustal  layers  and  the  associated  velocity  distributions 
as  reported  by  various  authors.  The  reason  for  this  is  probably  a  combina¬ 
tion  of  real  crustal  variations,  relatively  poor  sampling  densities  and 
the  non-uniqueness  inherent  in  the  interpretation  of  refraction  sei  mic 
data.  Consequently,  as  a  Cambridge  group  was  carrying  out  a  refraction 
experiment  across  the  North  Sea  in  the  summer  of  1980,  we  took  the  oppor¬ 
tunity  to  expand  the  shot-firing  scheme  in  order  to  obtain  a  densely 
sampled  land  profile  (CANOBE)  in  southern  Norway.  The  profile  was  posi¬ 
tioned  in  a  north-easterly  direction,  away  from  the  line  of  shots  off 
the  southern  coast  and  running  parallel  to  the  western  border  of  the 
Oslo  Graben  as  shown  in  Fig.  VI. 6.1. 

In  this  section  we  present  an  interpretation  of  the  high  quality  records 
obtained  along  the  densely  sampled  CANOBE  profile.  Synthesis  of  P-wave 
amplitudes  and  travel  times,  assuming  a  laterally  varying  structure, 
constitutes  the  basic  tool  for  mapping  the  crustal  structure  and  thickness 
in  the  coastal  areas  of  southern  Norway  and  the  general  features  of  the 
Moho  near  and  in  the  Oslo  Graben. 

Fie  1(1  work  and  prepa ra ti on  of  the  data 

The  CANOBE  proiect,  its  name  derived  from  the  participating  institutions 
Cambridge  University,  NORSAR  and  Bergen  University,  took  place  between 
the  26th  of  July  and  the  4th  of  August  1980.  A  total  of  11  recording 
instruments  were  available.  Four  explosions  from  a  coinciding  Cambridge 
North  Sea  refraction  proiect  (N2-N5  in  Fig.  VI. 6.1)  were  located  to  the 
east  of  the  Central  Graben  in  the  North  Sea,  and  were  used  for  the  land¬ 
recording  in  Norway.  In  addition  the  Royal  Norwegian  Navy  provided 
several  tons  of  AMATOL  in  the  form  of  torpedo  warheads  and  a  coast  guard 
ship  to  fire  shots  H1-H6  also  marked  in  Fig.  VI. 6.1. 


It 
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The  recurding  scheme  consisted  of  seven  legs  (Fig.  VI. 6.1),  one  per  shot , 
each  comprising  13  mobile  stations.  Legs  1  to  3  had  a  sensor  spacing  of 
4  km,  legs  4  to  6  a  6  km  spacing  and  1  eg  7  a  spacing  of  5  km.  beg  7 
transected  the  Oslo  Craben  just  north  of  Oslo  while  the  main  line  ran 
across  the  Precambrian  rocks  to  the  west,  extending  into  the  NORSAR  array 
siting  area.  Five  NORSAR  subarrays  (01A,  02C,  03C,  04C,  06C)  recorded 
the  shots  continuously,  thus  enabling  the  main  line  distance  range  to  be 
extended  to  about  315  km. 

The  initial  record  sections  were  filtered  between  0.2  and  15  Hz.  Fig.  VI. 6. 2 
shows  the  record  section  corresponding  to  the  main  line  for  shots  N5  and 
H1-H5,  reduced  by  8  km/s.  The  NORSAR  array  records,  from  450  km  onwards, 
have  been  low  pass  filtered  at  4.75  Hz  sampling  rate  20  Hz.  True  ampli¬ 
tudes  are  multiplied  by  distance  in  the  record  sections.  The  spectra 
for  typical  records  in  this  section  are  found  to  have  dominant  frequencies 
within  2.5  to  3.5  Hz. 

Interpretation  of  the  CANO BE  data 

Although  the  quality  of  the  data  and  the  relatively  dense  sensor  spacing 
along  the  profile  provide  and  excellent  basis  for  interpretation  we  are 
somewhat  encumbered  by  lack  of  observations  for  the  first  70  km  of  the 
profile  and  by  the  absence  of  a  reversed  coverage.  This  being  said,  the 
i nterpretational  results  are  as  follows: 

Interpretation  qf_  the  main  line  section  assuming  lateral  homogeneity 
Starting  with  Pg  and  Pn  velocities  taken  from  the  record  section  (Fig.  VI. 6. 2) 
and  Moho  depths  in  the  expected  range  of  30-35  km,  theoretical  time-distance 
curves  were  computed  until  a  satisfactory  fit  with  the  data  was  achieved. 

A  prominent  feature  in  the  record  section  are  the  large  amplitudes  of 
secondary  arrivals  at  about  186  km,  apparently  part  of  the  PyjP  branch. 

While  modelling  the  data  on  the  basis  of  travel  times  a  velocity-depth 
configuration  was  established  which  produces  a  focusing  effect  near  the 
outer  cusp  in  the  triplication  (Fig.  VI. 6. 3).  The  outer  cusp  is  made'  to 
terminate  at  250  km  hy  introducing  a  velocity  gradient  of  0.025  s“l 
in  the  depth  Interval  5.5  km  to  25  km,  and  an  increasing  velocity  gradient 


down  ti  32.5  km  where  a  velocity  of  8.1  km/s  Is  reached.  As  can  be  seen 
in  the  figure,  the  computed  subcritical  reflection  travel  times  do  not 
agree  with  the  secondar  '  arrivals  between  70  km  and  110  km. 

Now,  applying  the  reflectivity  method  to  the  model  in  Fig.  VI. 6. 3,  the 
above  gradient  is  sufficient  to  reproduce  large  amplitudes  in  the  retro¬ 
grade  travel-time  branch  of  the  observations  (Fig.  VI. 6. 4).  The  theoretical 
amplitudes  near  the  critical  distance  are  large  in  comparison  with  the 
observed  ones.  However,  further  attempts  to  model  these  amplitudes  in 
the  context  of  laterally  homogeneous  models  were  not  considered  due  to 
the  mentioned  travel  time  discrepancies. 

Interpretation  of  the  main  line  section  assuming  a  laterally  varying 
structure 

From  the  results  of  previous  profiling  experiments  and  seismological 
studies  in  southern  Norway  (Fig.  VI. 6.1)  crustal  structure  and  Moho 
depth  are  expected  to  vary  especially  around  the  endpoints  of  the  CANOBR 
profile.  It  is  essential  at  this  stage  to  deduce  a  model  for  the  first 
200  km  of  the  profile  for  which  PMP  observations  are  available  in  addition 
to  the  Pg  and  Pn  phases.  Close  examination  of  Fig.  VI. 6. 5  reveals  that 
strong  PjujP  amplitudes  are  confined  to  two  distinct  distance  intervals. 

The  first  is  around  115  km  and  is  considered  to  be  near  the  critical 
distance  where  relatively  large  amplitudes  are  expected.  The  second  is 
around  185  km  and  is  limited  to  a  few  seismograms.  The  latter  extreme 
P^P-amplitudes  can  be  attributed  to  focusing  effects  caused  by  waves 
touching  a  caustic.  While  modelling  these  large  amplitudes,  however,  the 
relatively  early  P^P  subcritical  arrivals  (Fig.  VI. 6. 3)  must  be  accounted 
for.  These  suggest  a  thinner  crust  as  compared  to  the  laterally  homogeneous 
model  of  the  previous  section  although  the  travel-times  of  the  Pn  arrivals 
from  150  km  onwards  must  still  be  retained. 

The  calculations  for  the  laterally  varying  models  are  performed  using 
Cassell's  (1982)  box  method  (based  on  zero-order  ray  theorv)  where  physical 
medium  parameters  are  defined  at  1  x  1  km  grid  points.  A  model  is  presented 
in  Fig.  VI. 6. 6  which  produces  good  agreement  between  the  theoretical  and 


observed  travel  times.  The  ma  i  n  features  in  the  model  consist  of  a  Molio 
with  .i  depth  ot  ’’ h  to  28  km  oft-eoast  which  increases  to  33  km  over  a 
distance  of  in  km  beneath  the  coastline.  The  velocitv  gradient  directly 
above  the  Moho  is  similar  to  that  in  Fig.  VT.8.3  except  that  it  does  not 
follow  the  topography  of  the  Moho  from  80  km  onwards  but  remains  at  a 
constant  depth  and  gradually  fades  away  after  200  km  as  shown  in  Fig. 

VI. h. 8.  In  the  coastal  area  the  main  discontinuity,  associated  with  a 
velocity  jump  from  8.8  to  7.8  km/s ,  occurs  at  a  depth  of  28-28  km  and  is 
overlying  a  velocity  gradient  zone  reaching  a  velocity  of  8.1  km/s 
at  34  km.  In  this  region  the  Moho  should  he  regarded  as  a  transition  zone 
in  which  its  depth  ’ s  not  clearly  defined.  Revond  110  km  the  Moho 
materializes  into  a  first  order  discontinuity  with  a  Pn  velocity  of  8.1 
km/s.  The  upper  mantle  velocity  gradient  is  approximately  0.008  s-^ . 

Synthetic  seLsmograms  were  calculated  for  Model  1  (Fig.  VI. 8. 7).  Note  that 
zero-order  ray  theory  as  used  here  is  not  exact  near  caustics  or  crictical 
points.  Nevertheless,  the  amplitude  distribution  in  the  synthetic  section 
qualitatively  justifies  the  choice  of  the  model.  The  amplitude  distri¬ 
bution  exhibits  distinct  peaks  corresponding  to  the  critical  point  of 
the  P^P  travel  time  branch  and  the  caustic.  The  first  theoretical 
amplitude  maximum  is  located  at  95  km  whereas  the  observed  maximum  is 
found  to  be  at  approximately  118  km.  This  discrepancy  is  not  surprising 
as  zero-order  ray  theory  does  not  include  wave  effects  in  the  critical 
region  beyond  the  critical  point.  A  similar  phenomenon  applies  to  the 
amplitudes  near  the  caustic,  which  are  exaggerated  by  zero-order  ray 
theory. 

The  Pn  arrivals  begin  to  undulate  at  distances  greater  than  270  km. 

The  interpretation  of  structures  in  this  distance  range  is  ambiguous  as 
we  cannot  determine  whether  the  travel  time  perturbations  are  caused  by 
Moho  topography  or  by  lateral  crustal  variations  or  a  combination  of 
both. 


In  Fig.  VI. 8. 8  , wo  possible  structural  configurations  are  presented  which 
satisfy  the  Pn  travel  time  data.  The  modelling  was  done  with  the  criterion 
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of  fitting  computed  travel  time  curves  to  the  trend  of  the  two  apparent 
velocities  in  the  Pn  arrivals.  These  are  a  model  with  a  varying  Moho 
(Model  1)  and  one  with  a  varying  surface  layer  thickness  (Model  2). 

Significant  constraints  are  imposed  on  the  interpretation  when  attempting 
to  account  for  the  late  arrivals  which  occur  between  16  and  20  s  in  the 
last  150  km  of  the  main  line.  These  arrivals  appear  to  be  surface  multiples 
corresponding  to  the  Pg ,  Pj^P  and  Pn  phases  of  the  first  arrivals.  When 
calculating  multiple  ray  paths  for  Model  1,  we  find  that  the  corresponding 
travel  times  fall  approximately  1  second  short  of  the  observed  arrivals. 
This  indicates  that  the  mean  crustal  velocity  is  too  high  as  these  waves 
travel  only  in  the  crust.  The  required  delay  is  achieved  for  Model  2  in 
Fig.  VI. 6. 8  by  introducing  a  Moho  depth  of  30  km  and  a  mean  crustal 
velocity  of  6.16  km/s. 

In  Fig.  VI. 6. 9  synthetic  seismograms  are  calculated  for  Model  2.  The  increa 
Ing  Pn  amplitudes  with  distance  are  in  general  agreement  with  those  in  the 
observations,  although  the  theoretical  ones  are  not  continuous  throughout 
all  the  records.  This  inconsistency  is  caused  by  unavoidable  lateral 
velocity  gradients  in  the  laterally  varying  model.  The  multiple  arrivals, 
however,  have  much  smaller  relative  amplitudes  than  those  in  the  observa¬ 
tions.  In  a  few  cases,  some  of  the  theoretical  records  lack  arrivals  due 
to  difficulties  in  finding  rays  to  each  receiver  location. 


Discussion  and  Conclusions 

In  this  study  we  have  demonstrated  that  the  use  of  amplitude  information 
as  a  supplement  to  travel  time  data  is  essential  in  delineating  eart  i 
structure.  The  outstanding  example  here  is  the  extremely  strong  amplitudes 
in  the  P^P  branch,  confined  to  a  relatively  small  distance  interval  around 
180  km,  which  appears  to  be  created  by  a  strong  velocity  gradient  in  the 
lower  crust.  In  this  respect,  the  calculations  of  amplitudes  for  2-0 
media  proved  very  important  as  compared  to  the  1-D  amplitude  calculations 
by  the  reflectivity  method. 
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The  Moho  is  found  to  he  27  to  28  km  deep  beneath  the  First  80  km  of  the 
profile  before  it  dips  downward  over  a  distance  of  30  km  beneath  the  coast 
and  subsequently  levels  out  at  a  depth  of  approximately  34  km.  The  inter¬ 
pretation  of  the  record  section  presented  is  only  approximate  for  the 
last  300  km  of  the  main  line.  Two  structural  models  were  proposed  for 
-this  distance  range.  In  order  to  keep  the  results  presented  here  consistent 
with  Kongsberg  seismograph  station  data  (Bungum  et  al  (1980)  find  a  Moho 
depth  of  34  km  at  Kongsberg  -  see  Fig.  VI. 6.1),  preference  should  be  given 
to  Model  l  in  Fig.  VI. 6. 8.  In  this  model  the  Moho  becomes  slightly  shallower 
as  it  approaches  the  margin  oi  the  Oslo  Graben  from  the  southwest.  After 
a  distance  of  310  km  along  the  profile  it  appears  to  sink  to  36  km  from 
where  it,  over  a  distance  of  180  km,  rises  to  a  depth  of  approximately 
35  km  beneath  the  N0RSAR  array.  This  Moho  depth  is  in  general  agreement 
with  Berteussen's  (1977)  Moho  depth  values  which  are  based  on  spectral 
ratio  analysis  of  long  period  P  waves  recorded  at  N0RSAR. 

On  the  other  hand.  Model  2  of  Fig.  VI. 6. 8  resulted  from  an  attempt  to  cor¬ 
relate  the  strong  secondary  arrivals  in  the  last  100  km  of  the  profile, 
interpreted  as  multiple  surface  reflections  of  the  ,  Py|P  and  Pn  phases. 

The  travel  time  delays  in  the  crust,  which  are  required  for  the  correlation 
of  these  multiple  phases,  infer  a  decrease  in  mean  crustal  velocity  by  as 
much  as  0.3  km/s  and  an  elevation  of  the  Moho  by  5  km  in  comparison  with 
Model  1.  Still,  the  synthetic  section  (Fig.  VI. 6. 9)  for  Model  2  was  not 
particularly  successful  in  matching  the  amplitudes  of  multiple  phases, 
and  the  possibility  remains  that  the  phases  in  question  have  not  been 
correctly  identified.  So,  in  order  to  conform  with  previous  investigations, 
we  adopt  Model  1  as  representative  for  the  main  profile. 

S.  Mykkeltveit 
F,  .S.  Husebye 

B.  Cassell  (Univ.  of  Cambridge,  UK) 
R.  Kanestrpm  (Univ.  of  Bergen) 
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Fig.  VI. 6.1  Simplified  geological  map  of  southern  Norway  with  CANOBE  shot 
points  and  recording  legs.  Results  from  previous  profiles  are 
reported  by  Sellevoll  &  Warrick  (1971)  and  Kanestr0m  &  Nedland 
(1975)  for  the  Flora-Asnes  and  Fed je-Grimstad  profiles,  Kanestrom 
&  Haugl and  (1971)  for  the  '3-4’-profile,  Tryti  &  Sellevoll 
(1977)  for  profiles  in  the  Oslo  Graben,  Weigel  et  a]  (1970)  for 
the  Skagerrak  profile  and  Mykkeltveit  (1980)  for  the  Arsund- 
Otta  profile.  Also  shown  is  the  NORSAR  array  siting  area  with 
the  original  22  subarrays,  each  comprising  6  short  period  in¬ 
struments.  Five  of  the  subarrays  (bold  rings)  recorded  the 
CANOBE  shots. 
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Fig.  VI. 6. 2  Main  line  record  section.  Amplitudes  are  multiplied  by  distance 
and  the  seismograms  are  filtered  between  0.2  and  15  Hz.  The 
records  beyond  450  km  were  recorded  by  the  NORSAR  array. 

The  continuation  of  the  crustal  phase  to  the  origin  Is  not 
clear  although  an  intercept  time  of  1.1  s  indicates  a  travel 
time  delay  in  the  upper  crust. 


Fig.  VI. 6. 3  Travel  time-distance  curves  for  the  laterally  homogeneous 

velocity-depth  model  (inserted).  The  portion  of  the  velocity- 
depth  function  corresponding  to  unrecorded  arrivals  in  the 
first  70  km  of  the  record  section  is  dashed.  Crosses  indicate 
observed  arrivals.  Note  the  discrepancy  between  observed  and 
theoretical  subcritical  PmP  which  indicates  an  elevated  Moho 
beneath  the  first  part  of  the  profile. 
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Fig.  VI. 6. 4  Synthetic  seismograms  computed  by  the  reflectivity  method  for 

the  model  in  Fi^.  VI. 6. 3.  The  source  signal  has  a  dominant  i 

frequency  of  3  Hz  and  4  extrema.  The  dot  indicates  the  outer 
cusp  in  the  P^P-surface  layer  multiple  caused  by  the  artificial 
first  order  discontinuity.  The  strong  P^P  amplitudes  are 
produced  by  the  velocity  gradient  in  the  crust.  Amplitudes  are 
multiplied  by  distance. 


Fig.  VI. 6. 5  Arrival  picks  in  the  main  line  record  section.  Amplitude  maxima 
in  the  secondary  arrivals  occur  at  113  and  IRS  km. 
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Fig.  VI. 6. 6  Comparison  of  observed  travel  times  and  theoretical  ones  computed 
for  Model  1  in  Fig.  VI. 6. 8.  The  vertical  radii  of  the  symbols 
indicating  the  observed  times  represent  the  reading  errors.  K 
denotes  the  Moho  depth  derived  from  Kongsberg  seismograph  station 
data.  (Bungum  et  al,  1980) 


Fig.  VI. 6. 7  Synthetic  seismograms  for  the  box  method  (Cassell,  1982).  The 

source  wavelet  has  a  dominant  frequency  of  3.1  Hz,  and  amplit  ides 
are  multiplied  by  distance.  Note  the  large  amplitude  at  176  kn 
resulting  from  the  velocity  gradient  in  the  lower  crust. 
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VI. 7_  Fennoscandlan  noise  survey 

During  the  summer  and  fall  of  1981,  the  first  part  of  a  noise  study  was 
performed  involving  various  sites  in  Finland,  in  northern  Norway  (Finmark), 
and  in  southeastern  Norway.  The  purpose  of  this  survey,  which  still 
continues,  is  to  measure  the  ambient  noise  level  as  a  function  of  frequency 
for  various  sites  in  Fennoscandia ,  and  to  measure  noise  correlation  as 
a  function  of  short  interstation  distances  for  some  of  these  sites.  Special 
emphasis  has  been  given  to  high  frequencies,  in  order  to  obtain  information 
relevant  to  the  construction  and  possible  deployment  of  small  aperture 
seismic  arrays  aimed  at  detection  and  analysis  of  regional  seismi'.  signals. 


Recording  equipment  and  logistics. 

The  main  field  equipment  used  has  been  Kinemetrics  PDR-2  'Compuseis'  recorders 
combined  with  Geotech  S-13  seismometers,  and  preamplifiers  were  needed  in  order 
to  obtain  sufficient  amplitude  resolution  for  the  background  noise.  The 
data  were  recorded  on  cassettes  and  subsequently  played  bad  using  a  Kine¬ 
metrics  CCS-1  playback  system  and  an  IBM  4341  computer.  A  r uch  faster  play¬ 
back  procedure  was  later  developed  based  on  a  PDP  11/34  computer,  end 
the  development  of  these  programs  was  an  effort  that  was  completed 
late  spring  1982. 

In  addition  to  the  data  recorded  by  the  PDR-2  field  recorders,  the  noise 
study  has  also  included  some  40  Hz  NORESS  data  recorded  at  the  NORSAR  data 
center.  These  data,  from  8  of  the  NORESS  seismometers,  have  been  tr  nsmitted 
in  an, '■log  form  to  Kjeller  and  recorded  there  via  the  now  MODCOMP/  IB!'  4331 
connection. 

Methods  of  analysis 

The  procedure  for  computing  power  spectra  depends  on  the  nature  of  the  data 
under  analysis.  Basically,  these  would  fall  into  three  main  categories: 

a)  Transient  signals: 

The  Energy  Density  is  estimated  as 
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Px(f)  =  At2X(f)X(f)  (nm2/Hz2) 

N-l 

where  X(f)  =  \  x(n)exp(-2nifn/N) 
n=0 

At  =  sampling  interval  (sec) 

N  =  number  of  samples 


b)  Stationary  signals  (noise): 

The  Power  Density  is  estimated  as 


At  - 

Px(f)  =  — -  X(f )X(f)  (nm2/Hz) 

N 

or  with  block  averaging  (M  blocks) 

At  M  - 

Px(f)  =  _  l  X.(f)X.(f) 

N«M  j=l  J 


c)  Periodic  signals: 

The  Power  is  estimated  as 

Px(f)  =  —  X(f)X(f)  (nm2) 

N2 


It  is  the  definition  under  b)  which  has  been  used  in  the  present  case,  for 
analysis  of  stationary  noise,  while  a)  would  be  required  for  analysis  of 
earthquake  or  explosion  signals.  Needless  to  say,  th<  normalization  procedure 
is  important  when  faced  with  the  problem  of  computing  spectra  with  reference 
to  absolute  ground  motion. 


Noise  correlation  p  (zero  lag  only)  has  been  computed  using  the  same  block 
lengths  as  for  the  spectral  estimates  and  averaged  over  the  same  number  of 
blocks.  The  estimates  are  considered  only  as  functions  of  interstation  distances 
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(azimuth  not  used),  and  averages  have  been  computed  over  certain  distance 
intervals.  Since  correlations  are  bounded  by  unity,  they  are  not  normally 
distributed  and  standard  deviations  can  therefore  not  be  computed  the  con¬ 
ventional  way.  This  problem  was  solved  by  transforming  the  correlations 
into  a  new  variable: 

z  =  i  I  *,n(l+p)-£n(l-p)  ] 

which  will  be  almost  normally  distributed  as  shown  by  Fisher.  Confidence 
limits  are  then  taken  as  z  =  z+oz  and  z  =  z-cz,  and  the  corresponding 
levels  in  p  are  found  by  the  transformation: 

p  =  tanh(z) 

Re  suits 

The  sites  which  have  been  analyzed  so  far  are  shown  in  Fig.  VI. 7.1  and  their 
names  and  locations  are  listed  in  Table  VI. 7.1,  where  also  the  main  results 
in  terms  of  spectral  levels  are  given.  For  most  of  the  locations  at  least 
2  time  intervals  are  presented,  with  results  given  for  7  frequencies,  i  ach 
separated  by  one  octave. 

It  should  be  noted  here  that  for  all  of  the  locations  from  1  to  7,  where  the 
PDR-2  field  equipment  was  used,  there  was  a  considerable  problem  in  finding 
time  intervals  that  were  not  contaminated  by  high  frequency  noise  of  more 
transient  nature.  The  sources  here  were  mostly  of  cultural  origin,  but  also 
In  some  cases  wind  and  effects  due  to  poor  coupling  between  seismometer  and 
ground  (for  many  of  the  sites  the  seismometers  were  placed  on  surface  exposure 
of  rock).  This  problem  is  of  course  connected  to  the  fact  that  we  are  inter¬ 
ested  in  very  low  levels  of  the  the  ambient  noise  (down  to  0.01  nm)  and  at 
frequencies  much  higher  than  with  more  conventional  seismometers  (sampling 
rate  62.5  Hz,  high-cut  filter  at  25  Hz).  For  the  future  measurements  we  will 
counter  this  problem  by  more  careful  siting,  concrete  pads  for  the  seismometers, 
wind  protection,  etc. 
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la  addition  to  the  first  sites  listed  in  Table  VI. 7.1,  where  the  PDR-2  equip¬ 
ment  was  used,  we  have  also  analyzed  some  intervals  of  40  Hz  NORESS  data 
as  mentioned  above.  The  reason  why  the  results  there  are  given  only  up  to 
4  Hz  is  that  the  linearly  ranged  quantum  unit  for  that  temporary  system  is 
too  large  to  give  sufficient  amplitude  resolution  for  higher  frequencies. 

Among  the  conclusions  that  can  he  drawn  from  the  spectral  results  are  the 
following : 

1)  The  spectral  level  at  1  Hz  is  for  most  sites  in  the  range  between 
0  and  10  dB  relative  to  1  nm^/Hz.  Although  this  is  10-20  dB  higher 
than  the  SRO/ASRO  low  noise  points  (Peterson,  1980)  and  the  Queen 
Creek  level  (Fix,  1972),  it  is  still  a  reasonably  good  level  as 
compared  to  most  seismic  stations  of  the  world. 

2)  The  slope  of  the  spectrum  falls  off  with  about  20  dB/octave  from  0.2b  Hz 
and  up  to  l. 5-2.0  Hz;  above  that  frequency  the  slope  is  around  10  dB/ 
octave  all  the  way  up  to  16  Hz  (see  also  Fig.  VI. 7. 2).  This  is  contrary 

to  the  stronger  reduction  in  slope  which  is  commonly  observed  for  frequencies 
above  about  5  Hz.  We  have  confidence  in  our  results  here,  however, 
since  i)  our  use  of  preamplifier  asures  a  good  quantification  for  the 
noise,  and  ii)  the  high  cut-off  filter  and  sampling  frequency  gives 
a  good  relative  response  at  these  higher  frequencies  (see  Fig.  VI. 7. 3). 

3)  The  noise  level  at  0.25  Hz  (which  is  the  location  of  the  well-known  micro- 
seismic  peak)  shows  much  stronger  variation  in  time  than  in  space.  It  is 
clear  that  the  noise  level  up  to  at  least  1  Hz  correlates  (although 

with  much  smaller  fluctuations)  with  these  4-second  microseisms,  which 
are  mostly  attributed  to  atmospheric  disturbances  over  oceans.  This  is 
consistent  with  the  change  in  spectral  slope  around  2  Hz  which  could 
indicate  two  separate  noise  processes  (Adair,  1982),  and  it  is  also 
consistent  with  the  results  of  Ringdal  and  Bungum  (1977)  who  studied 
long-term  noise  variation  within  different  frequency  bands  at  N0R8AR. 


J 
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4)  Although  there  ts  a  tendency  (as  expected)  that  the  Finnish  sites 
get  quieter  as  one  moves  eastward,  the  noise  levels  there  are  not 
significantly  better  than  those  observed  In  southeastern  Norway. 

The  northern  Norway  sites,  on  the  other  hand,  seem  to  be  slightly 
higher  in  noise  level,  which  should  be  expected  from  the  proximity 
of  these  sites  to  the  offshore  noise-generating  areas. 

5)  For  most  of  the  noise  data  analyzed  so  far,  there  is  no  significant 
difference  between  vertical  and  horizontal  components.  In  a  few  cases, 
however,  the  horizontal  components  chow  a  slightly  higher  noise  level. 

Results  in  terms  of  noise  correlation  are  presented  in  Figs.  VI. 7. 4-6, 
where  Fig.  VI. 7. 4  is  from  location  7,  covering  distances  from  50  to  600  m. 

The  scatter  in  the  results  is  fairly  large,  which  in  part  is  due  to  the 
fact  that  the  data  are  taken  from  non-simultaneous  measurements  (3  seis¬ 
mometers  were  successively  moved  around  to  12  sites),  and  in  part  caused 
by  the  previously  mentioned  high-frequency  transient  noise  disturbances. 

There  are  clear  indications  in  Fig.  VI. 7. 4  that  negative  correlation  occurs 
at  certain  distances  (dependent  on  frequency  bands),  which  was  shown  by 
Mykkeltv  !it  et  al  (1982)  to  be  consistent  with  a  model  of  propagating  no'se 
under  is  >tropic  conditions. 

Correlation  results  from  the  NORESS  array  sampled  at  40  Hz  are  shown  in 
Figs.  VI. 7. 5-6,  where  the  confidence  limits  now  are  much  smaller.  The 
tendency  for  the  noise  to  correlate  negatively  is  now  even  clearer,  although 
the  negative  deflection  is  not  as  pronounced  as  in  the  20  Hz  NORESS  results 
presented  by  Mykkeltveit  et  al  (1982). 

The  examples  in  Figs.  VI. 7. 5-6  are  picked  out  so  as  to  show  a  large  difference 
in  correlation  levels  in  the  first  frequency  band  (1-3  Hz).  It  is  interesting 
to  note  that  although  there  may  be  a  small  difference  also  in  the  2-3  Hz 
band,  there  is  no  difference  for  higher  frequencies.  It  seems  therefore  that 
the  main  correlation  variations  tend  to  occur  within  the  frequency  band 
covered  by  the  steepest  part  of  the  spectrum  (up  to  1. 5-2.0  Hz),  which  again 
point  towards  different  noise  processes  below  and  above  this  frequency. 


There  is  another  interesting  observation  in  the  fact  that  the  data  used  in 
Fig.  VI. 7. 5  (day  80/1982)  and  in  Fig.  VI. 7. 6  (day  130/1982)  have  reasonably 
similar  spectra  (see  Table  VI. 7.1),  in  spite  of  very  different  correlation 
levels.  There  is,  however,  a  difference  in  the  sense  that  the  spectrum  for 
day  80  (Fig.  VI. 7. 5)  peaks  at  a  higher  frequency  than  for  day  130  (Fig. 

VI. 7. 6).  This  points  towards  a  smaller  distance  to  the  main  noise  generation 
area  and  possibly  shallower  water  if  ocean-generated  (cf.  Bungum  et  al , 

1971).  In  fact,  the  1-3  Hz  noise  correlation  on  day  80  is  higher  than  for 
any  of  the  other  NORESS  intervals  included  in  Table  VI. 7.1,  which  all 
fall  between  the  levels  of  days  80  and  130.  It  should  be  mentioned  here  that 
some  correlation  measurements  obtained  at  location  4  are  consistent  with 
the  results  presented  above. 

The  main  conclusion  in  terms  of  noise  correlation  is  therefore  that  there 
is  a  fairly  good  stability  in  time  and  space  for  frequencies  above  about  2  Hz, 
while  there  is  some  variation  in  time  for  lower  frequencies.  This  variation 
correlates  more  with  the  shape  of  the  noise  spectrum  than  with  the  absolute 
level,  pointing  towards  differences  in  the  noise-generating  areas  as  the 
main  factor  of  influence  with  respect  to  the  correlation  properties  of  the 
noise. 


H .  Bungum 
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NQBSAP  NOISE  STuOy  SITES 


Fig.  VI. 7.1  Sites  for  the  Fennoscandian  noise  survey  1981/82. 
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Table  VI. 7. I  Name  and  location  for  the  8  noise  survey  sites  shown  in  Fig.  VI. 7.1. 

The  power  spectral  level  for  each  of  the  sites  are  given  for  7 
frequencies,  each  separated  by  one  decade. 
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Fig.  VI. 7.3  Relative  response  functions  for 

1.  NORSAR,  20  Hz  sampling,  5  Hz  filter 

2.  NORSAR,  20  Hz  sampling,  8  Hz  filter 

3.  NORSAR,  40  Hz  sampling,  12.5  Hz  filter 

4.  PDR-2,  82. 5  Hz  sampling,  25  Hz  filter 
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F1„.  vi. 7. 4  Noise  correlation  vs.  distance  for  location  no.  7  (see  Table  VI. 7.1) 
where  3  seismometers  were  successively  moved  around  to  cover  inter 
station  distances  between  50  and  600  m.  The  four  graphs  cover  the 
frequency  hands  1-3,  2-4,  3-5  and  4-6  Hz,  respectively. 
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Fig.  VI. 7.  5  Noise  correlation  vs.  distance  for  location  no.  8  (NORF.SS)  and 
day  80/1982  (see  Table  VI. 7.1).  Fight  channels  sampled  at  40  H 
have  been  used,  and  the  frequency  hands  for  the  four  graphs 
are  the  same  as  for  Fig.  VI. 7. 4. 
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